JOINT MOMENTS OF DERIVATIVES 
OF CHARACTERISTIC POLYNOMIALS 



PAUL-OLIVIER DEHAYE 
Pour Annie & Jean-Paul 

Abstract. We investigate the joint moments of the 2fc-th power of the char- 
acteristic polynomial of random unitary matrices with the 2/i-th power of the 
derivative of this same polynomial. We prove that for a fixed h, the moments 
are given by rational functions of k, up to a well-known factor that already 
arises when h = 0. 

We fully describe the denominator in those rational functions (this had 
already been done by Hughes experimentally), and define the numerators 
through various formulas, mostly sums over partitions. 

We also use this to formulate conjectures on joint moments of the zeta 
function and its derivatives, or even the same questions for the Hardy function, 
if we use a "real" version of characteristic polynomials. 

Our methods should easily be applicable to other similar problems, for 
instance with higher derivatives of characteristic polynomials. 

More data is available online, either on the author's web site or attached 
to the I^TjTjX source of this arXiv submission. 
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1. Introduction 

In Section 1.1, we merely define what is meant by joint moments of character- 
istic polynomials and state the results obtained in this paper. In Section 1.2, we 
motivate these Random Matrix Theory results by Number Theory questions and 
explain the interest of joint moments in the context of the Riemann ^-function. In 
Section 1.3, we discuss our techniques, which are essentially in Representation The- 
ory and Algebraic Combinatorics. The organization of this paper is summarized in 
Section 1.4. 



1.1. Presentation of results. We take for the characteristic polynomial of a x 
TV unitary matrix U 



N 



zu{o) ■■= n(i 

where the 9jS are the eigenangles of U. 
We define 



(1) 



Vu{0) 



It is easily checked that for real 9, Vu{0) is real and |V[/(0)| — \Zij{d)\. 

In this paper, we will investigate the averages (with respect to Haax measure) 



\M\j,{2k,r) 
iM)j,{2k,r) 

|VU(2fc,r) 

and their asymptotics 

|A^|(2fc,r) 
{M){2k,r) 
|V|(2fc,r) 



zum 



2k 



Zu{Q) 



U(7V) 



vum 



2k 



Vu{0) 



V{N) 



lim \M\j^{2k,r)/N 

N^OO 

lim {M)j^{2k,r)/N''' 



k^+r 



lim |VU(2fc,r)/iV'=' 



It is easy to show (by expanding the Haar measure explicitly) that the averages at 
finite N only make sense when 2k — r > —1. For the asymptotics, the normalization 
by N'^'^^ is due to Hughes [Hug05] (and proved in this paper anyways). 

This and related problems have been looked at by Conrey, Rubinstein and 
Snaith [CRS06], Hughes [HugOl, Hug05], Hughes, Keating and O'Connell [HKOOO], 
Forrester and Witte [FW06a], Mezzadri [Mez03]. However, much mystery remains, 
in particular for the dependency in r (when rgM \ N). 

While r e M \ N remains out of reach, we offer in this paper an alternative 
approach that uncovers some of the structure in those averages: 
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Theorem 1. For r G N and fc G C, the moments (Al)(2fc,r) are essentially given 
by rational functions, i.e. as meromorphic functions of k we have 

where Xr and are even monic polynomials with integer coefficients and with 
deg Xr = deg and G is the Barnes G-function [HKOOO, Appendix] . 
Moreover 

Yr{u)= W 

l<a<r-l 
a odd 

with the aa{-) given by 

—a + ^/ + Ar 
(^a{r) = 

We derive from this a similar result (Theorem 16, page 22) for \M\{2k,2h) and 
|V|(2fc, 2/i) (for h an integer). Finally, we have explicit expressions for {M.)(2k,r) 
given in Theorem 15, page 20 and Theorem 18, page 30 which allow us to com- 
pute the Xr{u)s, as given in Table 2, page 23, and additional data (available in 
Section 7). 

This structure had been guessed a few years ago by Chris Hughes based on 
computational evidence [Hug05]. The author is deeply thankful to him for freely 
sharing and explaining all of his previous unpublished work. 

This paper was initiated at Stanford University^ while the author was finishing 
his Ph.D., mostly worked on at Merton College, University of Oxford and finalized 
at Institut des Hautes Etudes Scientifiques, Bures-sur-Yvette. The author wishes 
to thank his hosting institutions for their support as well as his Ph.D. adviser, 
Dan Bump, and Persi Diaconis, Masatoshi Noumi and Peter Neumann for helpful 
discussions. 

1.2. Motivation. Ever since the works by Keating and Snaith [KSOOb, KSOOa], 
the Riemann ^-function can be (conjecturally but quantitatively) better understood 
through modelling by characteristic polynomials of unitary matrices. The classical 
example concerns moments. Let 



G(2fc+1) 



p prime m— 



p 



Then one can prove (fairly immediately, using the Selberg integral) that 
(3) |M|(2fc,0) = g(fc), 

which according to the Keating-Snaith philosophy leads to the following conjecture 
(for k > -1/2): 



(4) ^1 C(5+«) 







T 

dt g{k)a{k) ( lo; 



27r 



The main point is thus that a{k) is obtained by looking at primes, while g{k) is 
guessed at from the Random Matrix side. 

Observe also that Equations (3) and then (4) can be analytically continued in k. 



^with support from FRG DMS-0354662 
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Table 1. Summary of results on l{2k,2h) and J{2k,2h), when 
/i 7^ 0. The values for \M\{2k,2h) and \V\{2k,2h) are given as 
obtained from Theorem 16. Observe that the fifth column equals 
the product of the third and fourth. The "source" column refers 
to the paper where the result in the fifth column was published. 



k 


h 


a{k) 


\M\{2k,2h) 


X(2fc, 2h) 


source 


1 

2 
2 


1 
1 
2 


1 

6 

7r2 
6 

7r2 


{ 1\ 1'^ (2x1)^ 
V2V 1 (22-12)1 

/ 1 \ l2 (2x2)2 
V22 j 12 (42-12)1 

f 1\ l'^ (2x2)''-8(2x2)2-6 


1 
3 

2 

157r2 
61 


[Ing26] 
[Con88] 
[Con88] 


\2^ 1 12 (42-12)1(42-32)1 


1680^2 



k 




a(fc) 


|V|(2A:,2/i) 


J{2k,2h) 


source 


1 


1 


1 


/ 2! \ l2 1 
V1!2'V 1 (22-12)1 


1 

12 


[Ing26], 


see also [Hug05] 


2 


1 


6 


( 2\ \ 1^ 1 


1 


[Con88], 


see also [Hal02b] 


7r2 


\1\2^) 12 (42-12)1 


1207r2 


2 


2 


6 


4! ^ 1^ 1 


1 


[Con88], 


see also [Hal02b] 


7r2 


V2!2<i j 12 (42-12)1(42-32)1 


11207r2 



Many of the authors cited above have now shown that this philosophy should be 
extended to derivatives of characteristic polynomials. 

In particular, |A^|(2fc,r) should show up as the RMT factor of^ 



I(2fc,r) 



lim 



1 



,1 



2fc-r 



.1 



dt. 



3(^+it) 



2k- 



dt. 



and similarly |V|(2fc,r) is needed for 
J{2k,r):^ lim ^—r^ 

where 3 is Hardy's function (the relationship of 3 to C is analogous to the relation- 
ship of Vu to Zu, i.e. when t G E, 3(i + it) e M and ±3(i + it) = \C{1 + it)\ ). 
More precisely, it is expected that 

X(2fc,r) = a(fc) \M\{2k,r) and J{2k,r) = a{k) |V|(2fc,r). 

Thus, Theorems 1 and 16 give us a conjectural handle on the moments of C, and 3- 
One can compute some small cases (for integer k and r) and show that they 
agree with previous Number Theory (proved) results. This had already been done 
before and is repeated in Table 1. 

However, while Keating and Snaith obtained a full conjecture for T{2k, 0) and 
J'{2k,0) by computing |7W|(2fc, 0) and |V|(2fc,0), for the case of joint moments 
this goal remains elusive. All the available formulas for |A^|(2fc,r) or |V|(2fc,r) 
are rather inadequate. In particular, those formulas are limited to r :— 2h {h an 

■^It is a conjecture of Hall [Hal04] and Hughes [HugOl] that this is the appropriate normalization 
with respect to T. 
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integer), they are hard to compute for large values of k and h, they obscure some of 
the structure in the results, and finally they cannot be analytically continued in h. 

The analytic continuation would be important, because Conrey and Ghosh have 
proved in [CG89] that (on RH) 

^(2,1) = - 
and hence effectively conjectured'^ 

|V|(2,1) = 



47r 



An 

as well since a(l) = 1. In order to get this, we would need to have a sufficiently 
nice formula for \V\{2k,2h) that would allow for analytic continuation in h. We 
have simply been unable to do this but have no doubt that our results should be 
helpful for that goal (see the connection with Noumi's work below). 

On the other hand, the formulas obtained in Theorem 15, page 20 allow for much 
more effective computation than possible before, and we can compute longer tables 
for the different moments (see Section 7). 

This numerical data is useful as well, as Hall has devised (around 2002) a method 
that uses J{2k^ 2h) for all < /i < fc to produce a lower bound A(fc) on 

A := hmsup ^- , 



"^"^ log t„ 

where the t„ is the 7i-th positive real zero of C(l/2 + i^)- It is probably good to 
insist that this method does not depend on the Riemann Hypothesis, but only on 
values for moments! At the time of writing [Hal04], Hall only had the information 
he needed for A: up to 2 (conjecturally, up to 6). In Section 7, we present our 
conjectural data for J {2k, 2h) as a direct function of k for h up to 15 (see [Deh07] 
or the source of this arXiv submission for data up to ft, = 30). For a fixed h, various 
conjectural formulas are also given in this paper for J{2k, 2h) as a function of k. 
This, combined with Hall's method, should lead to more (conjectural) lower bounds 
on A. It is widely believed that A = oo so potentially we could also see if Hall's 
method has any hope to reach that, assuming only information on the J{2k,2h), 
but not on the Riemann Hypothesis. In other words, it would also inform us on 
the relationship between moment conjectures, the Riemann Hypothesis and the 
conjecture K — oo. We leave this to a further paper. 

Finally, Noumi in his book [Nou04] investigates the relationship between Painleve 
equations and expressions similar to one of the expressions we obtain for (A^)(2fc, r), 
in Theorem 18. Connections of this sort have been uncovered before (see [FW06a, 
FW06b] and works of Borodin), but an approach through Noumi's ideas would 
be original. One of our goals then would be to obtain analytic continuation for 
(A^)(2fc, r) in r, which would again allow to compute |V|(2, 1). We also leave this 
for further study. 

Our techniques are quite disconnected from the original motivation, so we discuss 
them separately. 

1.3. Techniques. As mentioned earlier, our techniques lie mostly in Representa- 
tion Theory and Algebraic Combinatorics. We look at the characteristic polynomi- 
als or the derivatives as symmetric functions of the eigenvalues of U , and express 
them in that way. We eventually express those symmetric functions in the most 
natural basis to use, the Schur functions. This basis is particularly suitable since 
those functions are also (irreducible) characters of unitary groups U(7V). We find 



•^This is completely backwards from the usual flow of conjectures from Random Matrix Theory 
to Number Theory, and possibly an unique instance of a reversal of this type. 
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ourselves integrating irreducible characters over their support (groups), which is 
very enviable! 

In order to express all the different functions in this basis of Schur functions, 
we use ideas present in a paper of Bump and Gamburd [BG06] and the author's 
thesis [Deh06b]. We will introduce those ideas as we need them. 

For a more thorough discussion of why a similar approach should always be 
attempted and other examples of its applications, please see the author's thesis 
and the results in [Deh06a]. 

Once we have a concise expression for the various moments, we still have to 
evaluate it. This will involve sums over partitions of values of Schur functions. 
After reparametrizing those sums over the Frobenius coordinates of the partitions, 
results of El-Samra and King were immediately useful to obtain the Schur values, 
and results of Borodin to handle the combinatorics of the sums. We then obtain 
a very big sum for the moments (Theorem 13), but that can directly be evaluated 
on computer (and thus checked against small N results) . After taking asymptotics, 
our results start simplifying into Theorem 15, enough to prove Theorem 1 on the 
general shape of those moments. However, the best expression is probably obtained 
once we use Macdonald's ninth variation of the Schur functions (Theorem 18). 

1.4. Organization of this paper. 

• In Section 2, we introduce all the non-standard notation we will be using. 

• In Section 3, we present the basic relations satisfied by the integrands 



Vu{0) 



1^^/(0) 1^'= 1 izum'" {WjY 1^^(0)1^'^ 

• In Section 4, we re-express the integrands as a sum in the Schur basis, in a 
way similar to Bump and Gamburd (via the Dual Cauchy Identity). 

• In Section 5, we engage in a long computation to evaluate the result ob- 
tained in the previous section, mostly using results of El-Samra and King, 
and Borodin. 

• Section 6 merely serves to tie what has been done in Sections 4 and 5 into 
the proof of Theorem 1 . 

• In Section 7 we present the data we are now able to compute, and partic- 
ularly discuss the position of the roots of |V|(2/c, 2h) in Section 7.2. 

• Section 8 describes two attempts to simplify our results further, one using 
Macdonald's ninth variation of the Schur functions, and the second imitat- 
ing a proof of the Cauchy identity. 

The bulk of this paper is contained in Sections 4 and 5. 



2. Notation 

We let N-i- be the set N\0. To avoid confusion with the index i, we have = — 1. 

We use V for a generic vector (of integers) (wi,--- ,Vd), and 'v for a sorted 
sequence of strictly decreasing integers vi > V2 > ■ ■ ■ > Vd, which we call a Frobenius 
sequence. Frobenius sequences are thus a special type of vectors. 

Sequences of weakly decreasing positive integers amount to partitions, and we 
stick with classical notation for those, i.e. A — (Ai, • • • ,Xi(x)), which defines 1{X). 
We also freely change our point of view to Young tableaux when discussing par- 
titions. Given a partition A of |A|, we denote its conjugate by A*. Define two 
sequences pi := Xi — i, qi := A* — i. They are strictly decreasing; Ai and A- are 
eventually 0, and hence pi — —i and qi = —i eventually. There exists d such 
that Pd > > Pd+i and > > qd+i- We call d the rank of A. We have that 
jj = (pi, • • • ,pd) and "q — (gi, • • • , qd) are Frobenius sequences, and we call jj 

and if the Frobenius coordinates of the partition A. We write A = ^ r- 
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Given p, we define CTp £ Sd such that sort(p) :— iPap{i)) is strictly decreasing 
(and hence a Frobenius sequence). This is thus not defined if pi ~ pj while i ^ j. 
We set sgn(p) := sgn((Tp), with the added convention that sgn(p) :— if Up is not 
defined. 

If A and p are partitions, A U /i is the partition obtained by taking the union of 
their parts. The partition {X^) has a.Y x X rectangle for Young tableau. 

We also use the notation [l^] for R copies of 1, used as argument to a (Schur) 
function. 



(5) 



3. Basic relations among the integrands 
We logarithmically differentiate Equation (1) to obtain 



Vu{0) 



and hence, when 6 is real. 



Zu{0) 



v{j{e) 



Vu{e) 
\Vu{e)) 

Zu{0) 



Zu{e) 



iV2 

2 

+ iN 



Zu{0) 



(6) 
(7) 
(8) 
(9) 



These basic relations give 
\M\^{2k,2h) 

\M\ {2k,2h) 

\V\^{2k,2h) 



\V\ {2k,2h) 



Y,(sNf-n .){M)N{2k,h+j), 

Y,x'^-^(%M){2k,h + j), 
i:(')(=f)""|A<U(2^.2i). 



These formulas are initially valid only when ft, is a non-negative integer, but the 
RHSs can be analytically continued by plugging in non-integer h and extending 
the sum to infinity^. We see thus that computing (A^)jv(2fc, r) would get us most 
of the way to \M.\N{2k,2h) or | V|Ar(2fc, 2ft,), and wc now focus on the integrand 



^Getting the correct analytic continuation can be tricky: The relation 

(10) |V|(2fc,2h) = . ) - (A^)(2fe,2h-j) 

j=o '-^ ' V^/ 

is also valid for integers h but here the RHS does not analytically continue in h to the LHS, since 
we exploit y^^^j = ( ) where it is critical that h be an integer. 
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4. Derivation into the Schur basis 

The goal here is to follow ideas similar to Bump and Gamburd's [BG06] in order 
to prove Proposition 4, page 9. One of their main tools was the Dual Cauchy 
identity. We encourage the reader to look at their first Proposition and Corollary 
for the unitary group, since this is all we really exploit from that paper. 

Lemma 2 (Dual Cauchy identity). If {xi} and {yj} are finite sets of variables, 

n 1 + x^yj = ^ S),t{xi)sx{yj), 

where the sum is over all partitions A and sx is the Schur polynomial. 

We (they) apply this Lemma setting {xj :— e'^^ : j G [1, • • • ,N]} to be the 
set of eigenvalues of U, and {yj := 1 : j G [1, • • • ,2k]}. We chose the notation 
sx{U) :=SA(e'^S--- ,e^^~). This gives 

J2^AU)sx{[^'']) = det(Id+C7)2'= 



,2k 



det{U) |det(Id+t/)|' 

,2k 



= S(fc«)(C/)|det(Id+C/)r 
or (replacing U by —U) 

\Zu{0)r = \det{Id~U)f' = (_i)'=^,^,„^([/)^(-i)|Ali^,,([l2'=]) 



We can also re-express 



Zu{0) 



i=i 

N 



J2 ' ■ 

j—1 m—1 



(11) = i lim_ J2 ^""PmiU) 



m—1 



where Pm{xi, - ■ ■ ,xn) is the m-th power sum + ■ ■ ■ + x"^ and we have used the 
same convention as for s\(U) of inputting the eigenvalues. We will use the same 
convention soon for the power sums px :— Yii P\i ■ 

In practice, we want the reader to just ignore the variable z and set it to 1. This 
will be justified a posteriori. 

Putting everything together, we thus get for |Z(7(0)|^'' ( z^(o) ) 
(12) {~\f^s^,.^{U)(xY.Prn{U)\ Y.i-ir^sxill"'])^^). 

\ m=l / A 

At this point, we will soon want to use the fact that the sx are characters of 
unitary groups. 

Indeed, li U £ U(A^) then when l{\) > N, we have^ sx{U) = 0, but when 
/(A), i(/x) < N, we have 

{sx{U)j;jU))^^^^-Sx,, 

i.e. for large enough N, sx is an irreducible character of U(A'^). This orthogonality 
is obviously good for our purposes, but the only obstacle is the need to express 



This is a consequence of the fact that sx(xi, ■ ■ ■ , Xn) = if 1{X) > n. 
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S{k") (U) (X]m=i Pm{U)Y exclusively in terms of Schur functions. This can be done 
and will require the Murnaghan-Nakayama rule. 

Let a ribbon be a connected Young skew-tableau not containing any 2 x 2-block. 
If a ribbon contains m blocks, it is called a m-ribbon. A first approximation to one 
version of the M-N rule says that s\pm is given by a signed sum of s^s, where ^ 
runs through all partitions obtained by adding a m-ribbon to A. 

If we average Expression (12) over U(A^), we could thus see A as running through 
all partitions obtained by adding r ribbons to the rectangle I^N'^') (this uses the fact 
that this lax version of the M-N rule is invariant under transpositions, since we have 
yet to discuss the signs). There are more conditions however. We also need /(A*) < 
N (since otherwise s\t{U) = 0, as in footnote 5), and we need Z(A) < 2k (since 
otherwise sa([1^'^]) — 0, again just as in footnote 5). In other words, A contains 
(iV'°) but is contained in l^N"^^). There are only finitely many (ways to obtain) such 
partitions, which will make the sum over As finite, and thus only finitely many sets 
of lengths of the r ribbons will contribute. This justifies a posteriori setting ^ to 1 
in Equation (11), but only when we can apply the dominated convergence theorem. 
This will only occur if we know of a bound on the integrand independent of z that 

whenever this is integrable, i.e. 



l2fe 



Z'uiO) 



Zu{0) 



is itself integrable. We can pick |Zt/(0)| 
only when 2k — r > —1. 

We now state a more precise version of the M-N rule 



Theorem 3 (Murnaghan-Nakayama). Let X be a partition and p be a vector with 
|A| = ^iPi- If is value of the irreducible character of S\\\ associated to A 
on the conjugacy class of cycle-type sort(p), then 

(13) Pp - 

A 



and (more importantly) 

(14) = E(-l)''^'^ 

s 

summed over all sequences of partitions S = (A'-'^-', A^^-*, • • • ,A^''-') such that r := 
1{X), = A(o) c A(i) C • • • C Af"") = A, and such that each A^^^ - A^^-^) is a ribbon 
of length p^, and ht{S) ^ Eiht(A(*) - A^^^^^). 

We have not defined the height ht of a ribbon, but rather than doing so or de- 
tailing the computation here, we only expose the idea: Equation (13) tells us that 
iJ2mPmY can be computed using the character values of symmetric groups, which 
can be evaluated by summing over sequences of partitions (A'"\ • • • , A*^*"'). For each 
such sequence, the sequence (A(°\ • • • , A^''^), with A^*^ := (iV'^) U A^'\ would be as- 
sociated with the combinatorics of the expansion of the product in Expression (12). 
Indeed the combinatorics of ribbon is unchanged under translations (down by k) as 
long as the partitions are kept within a rectangle (actually, a horizontally bounded 
region) . 

If the computation is explicitly carried out, we get the following result: 
Proposition 4. If 2k ~ r > —1, we have 

(15) (M)^(2fc,r) ^ i-iyJ2 E 4^(^^>ua([i'1), 

TIGN'^ a within 
+ kxN 

with the understanding that = i/ |A| 7^ fi^. 

For this result, we have preferred to index all the partitions containing (iV*^) but 
contained in (iV^fe^ ^ ^^fe^ y fo^. ^ (^j^ky 
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We are now left with the task of evaluating the RHS in Equation (15), which 
will turn out to be a tedious process. 

5. Main computation 

We are left with two problems. The first one is due to the characters of the 
symmetric group. Those are of course desperately hard to evaluate directly and 
individually. We are helped here because we will actually only evaluate something 
close to 

for given A. This amounts to computing the sum of values of the character over 
permutations with I cycles. The second issue is evaluating s^/^rfc^uA 

([l2fc]). The 

author had previously used the Weyl Dimension Formula to do this (see [Deh06b]). 
A formula giving that dimension in terms of the Frobenius coordinates of A is 
probably better adapted for our purposes. 

In addition, both "problems" combine extremely well, in that both expressions 
should involve a sign, which turns out to be the same. 

We will then sum our terms over all partitions, expressed in Frobenius coordi- 
nates. This amounts to summing over possible ranks (1 < d) and then pairs of 
Frobenius sequences of length d. 

5.1. The value of the Schur function in Frobenius coordinates. 

5.1.1. Dimension formula in Frobenius coordinates. El-Samra and King [ESK79] 
use the notation for s f-p"i ( [l''^] ) • 

l"q J 

Assume has d Frobenius coordinates. They prove that 



{R + Pi)\ 



{R-q,-l)W■q]^.{p^ + qJ + l) 



dxd 



i—l l<i<j<d — 1 



where the first expression is also known as the reduced determinantal form (cf. 
Foulkes [Fou51], as cited in [ESK79]). 

It is a consequence of Cauchy's Lemma that the two expressions in Formula (16) 
are equivalent: 

Lemma 5 (Cauchy). 
(17) ^ 



Pi + qj + 1 



= n (P' "Pj)(* - ii) n 



Observe that Formula (16) is positive (as it should, given that it is also a dimen- 
sion) because the pi and qi are strictly decreasing. 

However, the RHS of Formula (16) still makes sense if we plug in unsorted vectors 
p, q (with even the possibility of i ^ j but pi — pj). Hence this can be used to 

define sj-p-i ([l^]) as well, which is then skew-symmetric in both the piS and the 
IqJ 

QiS separately. This can be written 

(18) s.p,([l^]) = sgn(p)sgn(q)sr_,(p).([l«]). 

IqJ \soit(q)j 
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Observe that Formula (18) is still valid when sort(p) or sort(q) is not defined (this 

happens when two of the entries of p or q are equal) thanks to sgn(p) sgn(q) = 

(see conventions in Section 2)! 

Finally, it is helpful to remark that Formula (16) for s j-p-, ([l^]) can be seen as 

_ _ _^ _ _ _ _ _ 

a product indexed by the sets p U q and pairs in the set (p x p) U (q x q) U (p x q) . 

5.1.2. Evaluation of S(^iqk^tjx 

([P*^]). We take A = ||} to have d Frobenius coor- 
dinates. 

In total analogy with Equation (18), we first extend the definition of s^jv'')ua 
and set 

^ ' l_sort(q)J 

with the understanding (as before) that the value of the RHS should be if — pj 
(resp. Qi = Qj) for i ^ j. Again, this is skew-symmetric in the piS and separately in 
the QiS. 

We have the following Lemma 
Lemma 6. Let p, q be vectors with d coordinates. Then 

(19) V^)u{P} = 

J^l^fJ(p. + fc + i)W(7v + ,. + i)wj-^{i}(^^ 

Proof. By skew-symmetry, we really only have to check this for l-g'l' If we want 

to use Formula (16), we should look at the Frobenius coordinates of {^N'^'j U A. 
This would be rather unpleasant (particularly because the number of Frobenius 
coordinates would change for fixed N and k according to the A considered). 
Let us look instead at: 

^ -.^ {N + k~l,--- ,N), 
y := (2fc-l,... ,fc), 
c? := U I? (sorted) 
and (3 := U "q (sorted). 

Then (J and (3 are strictly decreasing, so those are Frobenius coordinates. The 
partition corresponding to those coordinates is obtained geometrically by sticking 
a (/c^'^) block to the left of (TV*^) U A, or equivalently to shifting (^N'') U A by fc spots 
to the right, while considering A = (Ai, • • • , A^) to have exactly k parts (with some 
possibly empty). 

Because of this, we have (as in [BG06, page 6]): 

Additionally, are the Frobenius coordinates of {{N + A:)'') U {k'^). Hence, 

for the same reason as above, we have: 

■^{1} ( [l^''] ) = ^HN+muik") ( [1^''] ) 
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When evaluating the product described in Equation (16) using the 1$ and /3 

^ — > ^ ^ — > — > 

coordinates, we have a big product taken over the sets a,/3,a x a, (3 x /3 and 

(J X /3 . We expand those index sets using = 1? U jj and /3 = U ^ . 

One can see that the products indexed byp,q,pxp,qxq and p x q 

together give S|^|([l^'']) = sx{[l'^^]). 

Similarly, the products indexed by 3c, y^, Icxlc, y^xj?' and x* x y^ give 

We are left with only "cross-products" to evaluate, for the index sets ^ x jj , 
Ic X "q , y^ X jj and y^ x "q . The definitions of 1? and ~y now give the result. □ 

5.2. Sums of characters over conjugacy classes with same number of cy- 
cles. Assume / (^|^|) is a function of pairs of vectors of the same length (say d). 

One can set /(A) := / ({|}), where A = {|}. 

The goal in this section is to evaluate sums of characters of the general form 

E ^^/(^)- 

We will eventually take /(A) = ^{n>')u\{\}-'^'^\) but there is no reason to limit 
ourselves in that way for a while. 

We rely on a few results of Borodin that give a slightly different version of the 
Murnaghan-Nakayama rule. 



5.2.1. Definitions. This is based on [BorOO, around page 15] and [Bor98, around 
page 6]. The relevant definitions (not included here) are fragment, the different 
block types, the filling numbers, filled structure, sign of a structure. 

Theorem 9 is almost in Borodin's work, and his definitions are used in Proposi- 
tion 10. Both of those results are used for Theorem 13, which can be read without 
looking at Borodin's papers. 

However, the first condition to have a fragment needs clarification in both papers: 
change 

(1) there is exactly one hook block that precedes the others 

to 

(1) there is exactly one hook block in each fragment. That hook 
block precedes any other block in the fragment 

We also would like to correct a statement in [BorOO] , in that linear horizontal or 
vertical blocks are positive, not just non-negative integers (in agreement with the 
other cited paper of Borodin [Bor98]). 

We can highlight one of the definitions: Any filled structure T with d fragments 
produces a set of pairs 

{(Pi,9i), • • ■ , {Pd.qd)} 

which consists of the filling p- and q- numbers of the fragments. 
The sign of T is defined as follows: 

sgn(T) = sgn(p)sgn(q)(-l)5:9-+-(^), 



where, as a reminder, the sgn inside the formula is if pi = pj (resp. qi — qj) for 
i 7^ j- 
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5.2.2. Simplified Mumaghan-Nakayama rule. Although we haven't defined any- 
thing, we state Proposition 4.3, taken from the first paper of Borodin: 

Proposition 7. For any two partitions A and p with |A| = \p\, we have 

T 

where the sum is taken over all filled structures of cardinality p — (pi, • • • , pi) such 
that the sequences (pi, • • • ,pd) and (qi, • • • , qd) of filling p-numbers and q-numbers 
of the structure T coincide, up to a permutation, with the Frobenius p-coordinates 
and q-coordinates of the partition X (i.e. X — |sort(^)} )■ 

The proof of this Proposition is quite simple: going back to the original pre- 
sentation of the Murnaghan-Nakayama rule in terms of hooks, Borodin analyzes 
what happens to Frobenius coordinates when subtracting hooks/ribbons. Each 
such subtraction corresponds to a block. There are three cases to distinguish: 
the hook/ribbon can be above or below the "Frobenius diagonal" or even overlap 
it. Those cases correspond respectively to linear horizontal blocks, linear vertical 
blocks, and hook blocks. 

This Proposition, as stated in Borodin's work, is slightly restrictive: there is no 
need for p to be a partition. Let p — (pi, pi) to he a, vector of positive integers 
and define (just as in Theorem 3) Xsort(p)- Then, by summing over all vectors 

P, we get 

Proposition 8. For any partition X, we have 

T 



+ 

where the sum is taken over all filled structures T of I blocks and with filling p- 
numbers (pi, • • • ,pd) and q-numbers {qi, - ■ ■ ,qd) such that X — |sort(^)}' 

Observe that d, the rank of A, has to be less or equal to I in order to have a 
structure. 

We now state the main Theorem we will use that is originated in Borodin's work. 

Theorem 9. Assume f is skew- symmetric within its two vector entries (sepa- 
rately), z.e. / ({-J^gj}) = sgn(p) sgn(q)/ ({|}) . Then, 

E E4/(A) - E E /({I}) Em)^^""^^"^ 

X within peW^ d=l pe[o,JV-l]'' T{p,q) 

^ q6[0,fc-l]'' 

where T(p, q) goes trough all filled structures of d fragments, I blocks, v{T) vertical 
blocks with filling p-numbers (pi, ■ ■ ■ ,pd) and q-numbers (<Zi, ■ • • , qd)- 

Proof. We start by summing Proposition 8 over As fitting inside a, k x N box: 

E E4/W = E E (-i)^''^"^'^^«g^(p)«g^(q)/(^) 



A within /)f=N' A within T(p,q) 

A.-xA^ ^ + kxN ^ ' 

= E 

A within T(p, a") 

kxN 



where the second sums in each RHS are taken over all filled structures T(p, q) of 
I blocks and d fragments such that the sequences of filling p-numbers (pi, • • • ,pd) 
and q-numbers (91, • • • , qd) of the structure coincide, up to two permutations, with 
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the sequences of Frobenius p-coordinates and g-coordinates of the partition A (i.e. 
A = Note that d changes with A. 

We then obtain the final result by seeing the double sum over A then permuted 
Frobenius coordinates of A as a sum over all vectors of appropriate lengths. 

We should not be concerned about vectors having two identical coordinates (say 
Pi = Pj), since the corresponding term in the RHS vanishes by skew-symmetry 
of/. □ 



5.2.3. Counting structures. We now need to compute the sum 

T(p,q) 

which is taken over the structures described above, i.e. for given I, d, p, q, v. 
It would help to know how many structures there are for each choices of those 
parameters. We prove the following Proposition. 

Proposition 10. There are exactly 



(20) #Til,d,p,q,v) = 



E 



d+J2si+ti=l 



{sd + td + --- + si + ti + d)\ 

n n ti\ ntiid + 1 - » + e + tj). 



n 



structures with d fragments, I blocks, filling numbers p = (pi,-'' ^Pd) and q = 
(9i) ■ ■ ■ ! Qd) o-nd V vertical blocks. The indices in the sum Sj (resp. ti) count hori- 
zontal (resp. vertical) blocks in the i*'^ fragment. 



Proof. This is a purely combinatorial problem. Given the number of vertical blocks 
on each fragment, we essentially have a partial order on blocks that we want to 
extend to form a linear order (across fragments). Part of the rules in the initial 
partial order say that the hook-block in the ith fragment precedes any other block 
in that fragment. We then need to fill the structure (i.e. choose filling numbers for 
each block). 

We can reverse this process: 

• We first choose the numbers of horizontal and vertical blocks and ti on the 
i*'^ fragment. We have the conditions that = ^(^) d+J2si + ti = I 
(i.e. there are I blocks in total, d hook, Si horizontal in the i^^ fragment 
and ti horizontal in the i**^ fragment). 

• Starting from the d^^ fragment, we decide where to insert the horizontal and 
vertical blocks of the i*^ fragment in the partial order that is established 
so far on the set of fragments from the i + 1st to the d^^ one. 

• We decide how to cut up the fragment into filled blocks, respecting the 
number of horizontal/vertical blocks decided upon earlier. 

The equality in the statement is intended to reflect clearly the layering described 
above: the sum corresponds to the first layer, while the other two layers correspond 
to one square-bracketed factor each. 

Observe that the relation Sd + td-\ \-si+ti + d= I could be used to simplify 

the numerator in this expression. 
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The only hard part is to derive for the second step 

isd + td + --- + si + ti + dy. _ 



{sd + td + --- + si + ti + dy. 

( (sd+td+i)x ^ 

■QgJJI^.ll (sd+*d+Sd-i+td-i+2)x 

\(sd+td+Sd-i+td-iH hsi+ti+d)y 



This is obtained by simphfying 

n 



i=0 



where the i*^ factor in the Jli^o"?^'^^^''^ counts the number of ways of choosing 
the hnear order on the blocks of the d — i"^ fragment, knowing the linear order 
restricted on the blocks of the fragments d — i + 1 to d. 

The first binomial factors intersperses the set of blocks of the d — j*^ fragment 
among the blocks of fragments d — i + 1 to d, while the second factor decides which 
blocks are horizontal and which are vertical. □ 

We wish to insist on the fact that the summand in Equation (20) is not symmetric 
in the p^s or the g^s, because the factor in the denominator niLi(^+l~*+X]j=i + 
tj) is not symmetric in the SjS or the tjS. For instance, Sd appears d times while 
s\ appears only once. 

5.2.4. Sum of determinants. We aim now to put together all the results obtained 
so far in this section, but we first need a quick lemma. 

Lemma 11. Let s and t be vectors of integers. Then, 

(21) E — ^^^^^^ = 



n (s^-sjKt^-ti) n 



1 + Si + t-j 

l<i<j<d l<ij<d ■' 

Proof. The proof proceeds as for the classical computation for the Vandermonde 
determinant: the LHS is skew-symmetric in s and t separately, has obvious poles 
as prescribed in the RHS (when Si^ + tjo — ^ 1 ) > ^^"^ the degrees in the RHS are 
appropriate. Up to a constant of proportionality, both sides are thus the same. 
This constant is shown to be 1 by looking at rates of decrease when si goes to 
infinity. □ 

Proposition 12. Assume f is skew- symmetric within its two vector entries (sep- 
arately), i.e. / ({^^Sg}) = sgn(p) sgn(q)/ ({p}) . Then, 



I 

E E 

TieN'^ X mthm d=l pg[o,iV-l]'' 

^ qe[0,fc-l]'' 



V t GN * \ / \ / J l<'i<j<d i<i,3<d ■> 
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Proof. We go head first and combine Equation (15) and Theorem 9: 



TtgN'i ^ within 

+ kxN 



qe[0,fc-l]'* 



E 



nfe:)(?:) 



^!E E 

<i=l pe[0,A'-l]'' 



E 



s.teN'' 

i:u=v 

d+J2s^+ti=l 



V 

sgn(p) sgn(q) Jlf C;) (^0 

n .'^.i n nti(rf + 1 - * + e + 



^!E E 

{1} 

within kxN 



+v 



E E 



It is now crucial to observe that for fixed p , q , s , t , the sign of the numerator of 
the summands (bracketed) will depend on the parity of a and t. Hence we obtain 



"t E /({|})(-i)^'-"- 



within fcxAT 



J2 E sgn(^)sgn(0)n 



d 



d+Y,s^+ti=l 



sgn(cr) sgn(T) 



V 
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=1 



{1} 



s , t eN° 

within fcxAf d+X]si+ti=' 



E 



sgn(a)sgn(r)n^ Cr;,)!^) 



The last hne is now perfectly cut for substitution using Lemma 11. After changing 
the range of summation on within fc x to p G [0, iV — l]'^,cl G [0,k — l]"^, we 
obtain the announced result. 

Admittedly, this is not very enlightening. It is thus worth highlighting what 
happens: the sums we are dealing with initially are sums over partitions. By using 
Frobenius coordinates, and sorting the partitions by their rank d, we are expressing 
the main sum into a sum over d of multisums in d variables. We thus now have sums 
over two sets of d strictly decreasing variables (the sets p and q ) of different ways 
of building up this partition (the data encoded in s and t). Using skew-symmetry, 

we can unsort the variables p and q to p and q and decide instead to sort the 

_ _ ^ — > 

variables according to "building blocks" , i.e. switch from s and t to s and t . □ 

5.3. Putting everything together. We combine all the information obtained so 
far, and simultaneously clear the restriction d + ^ + t.; = I in Formula (20) by 
encoding all the moments at once into an exponential generating function. 

Theorem 13. For a fixed k G N. the two series 



OO 



(22) J2(M)N(.2k,r)^-l- arid .(Ar.)([l'1)E E 



r>0 



S,\tj\{l + S, + tj) 



E 

dxd — 



qe [o,fc-i]'' 



{N + qj + l)(k){l+p^ + qj) 



dxd 



have equal coefficients of for r <2k + 1. 

Proof. A first necessary remark is that as a formal power series, the second series is 
well-defined: the sum to obtain the r^^ coefficient in that series reduces to a finite 
sum (because Si < pi and tj < q.j). 

We know from Equation (18) that s j-p-, ([l^*^]) is skew-symmetric in p and 

q (separately). Hence we can combine Equations (15), (16) and (19) and Proposi- 
tion 12 to obtain a huge sum. The main statement then follows from recombinations 
of the main product into determinants, using Cauchy's Lemma 5. □ 

Remarks on Theorem 13 

• This is a hypergeometric multisum (at least for fixed d), when we expand 
the determinants using Cauchy's Lemma. However, not even small ds seem 
tractable on computer. 

• A definite advantage of this formula is that it can be tested at finite N (by 
expanding the integral defining {Mff)(2k,r) symbolically using the Haar 
measure). This is helpful to confirm the results obtained so far. 

• We wish to insist on the idea behind this Theorem: initially we had a com- 
binatorial problem on structures (see Formula (20) that had no symmetry 
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for its summands in the s^s (resp. tis). We have exploited some skew- 
symmetry in the as and 6s in Formula (19) to change this. In particular, 
we have now switched from a sum over jj , "q, s, t to a sum over p, q, ~s*, t . 
We have also simplified the denominator in Formula (20). 
As a consequence of the previous point, we can now assume that the s^s 
are all different. The same is true for the tiS. 

This has useful consequences, especially for computational purposes. It is 
interesting to compute a bound on r such that partitions with d fragments 
will have a non-zero contribution to the final sum in (A4)jv(2fc, r) . We have 
r > d + '^Si + ti, and the s^s (resp. tis) should be all different. We can 
take them to be 0, 1, • • • ,d-l. We thus have r>d + 2^^^^ = d^. 



We now define 
(23) iJ^'*^-'^'* 



II ^ fc(iV-p) ('=)(- 1)9 fk + p\fk-l\fp\fq 



ge[o,fe-i] 

where the RHS is taken to be similar to the entries in one of the determinants in 
Equation (22). 

I have not been able to obtain a much better expression for this with Mathemat- 
ica. Normally, the package MultiSum [Weg] should be able to deal with multiple 
hypergeometric series, but this particular one is too complicated. We will thus fo- 
cus on an easier problem from now on, the problem of asymptotics (i.e. we switch 
from (Al)Ar(2fc,r) to (Al)(2fc,r)). 

5.4. Asymptotics. We need to compute asymptotics for }J^'^'^'* more precisely. 
Proposition 14. For a fixed integer k>\, when k > t, 

TTN.k,S,t 

(24) H''-''' := lim 

(25) = k J2 



Tik + i)T{s + i + t+l) 
T{i + l)T{k + s + t + i + 2) 



(26) 
(27) 



i=0 

^U=k 



1 T{2k ~ t)T{k + s + I) 
1 + s + tT{k - t)T{2k + s + 1)' 



fc(jV_p)fc(_l)g //fc-l\ /q\ 
iN + q+l)^l+p + q) k\ \ q )^ \t)' 

ge[0,fe-l] 

i.e. i^^'*^''*'* stripped of some of its terms of obviously lower order in p, N and q 
combined. We do this because we want to compute the leading order of H^'^^'^-'^ and 
there will be lots of cancellation due to the sum over q (as showed by Equation (28)). 
We thus wish to compute hmjv^oo #^.'=.«.*/7Vi+^+* = limAr^oo H^''''''* /N^+'+*. 



Proof. Define 

jjN,k,s,t _ ^, J2 



'Observe that Expression (27) is well defined, thanks to the bound k > t. 
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The proof of the equahty (24)- (25) essentially follows from two basic identities 
on formal series: 



{l-rX + r^X^ -...Y{l-sX + s^X^ 



^(-i)^5:(^+;-i).v-^x^ 



and 



i<k-i V y / 



(28) >' (-!)'( ^]q'X^ = {-lf+^{k-l)\X^-^ 



0<j 
0<q<k-l 

We expand the definition of iJ^'^-''^-* as a power series in q. The first identity 
indicates that we should only look at the coefficient of q''^^, which we obtain by 
using the second identity (set r := 1/N, s := l/{p+ 1)). We then let N tend to 
infinity, so the sum over p becomes a Riemann sum. Its limit is a /^-integral, and 
thus a /3-function appears, which can be expanded into a product of F-functions, 
giving the first equality. 

The equality (26)-(27) is immediate and is the only one to require the bound 
k>t. 

For equality (25)-(26)^, we first define 

T{k + a + i)T{s + a + t + t + l) 



T(a + i+ l)T(k + s + t + a + i + 2) 

2—0 



kT{a+k)T{a + s + t+l) f i,a+k,a+s+t+i 

T{a + l)F(a + fc + s + t + 2) ' ' [a+i,a+k+s+t+2^ ^ j . 



which satisfies H'^'^'^ = —H^,'^l . The second equality is merely a consequence 
of the definition of 3F2. 
Since (see [Mat]) 

/ X c + d-e + l , F(c+rf-e + l)F(e) 



(c - e+ l)(d - e + 1) V r(c)F(d) 



we get 



^k.s.t ^ 1 aT{a+k)T{a + s + t+l) 



l + s + t\ F(a + l)F(a + k + s + t + l)_ 
which lets us prove equality (25)-(26) using the relation H^^'^* = H^^'^' - H;-ff . □ 

Let G(-) be the Barnes G-function [HKOOO, Appendix]. It is a quick consequence 
of the Weyl dimension formula (see [BG06, Equation (18)] that 

^ n.2ki. G{k+lf ,2 
We use the previous Proposition to give a relatively concise expression for {A4){2k, r). 



This equality was first proved using Mathematica. Paul Abbott observed that the hyper- 
geometric function that appears is Saalschiitzian and extracted the following proof by tracing 
Mathematica's output. 
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Theorem 15. For a fixed fc e N, the two series 



(29) Y.^M){2k,r) 



{izY 



r>0 



and 



Gik + 1) 



(30) 



G(2fc+ 1) 



G(fc + 1)^ 
G{2k + l) 



1 




jjk,Si.tj 


sjtj!(l + Si - 


M,) 


dxd 




OO 

-E E 




1 






+ 


Sj + ij) 



dxd 



dxd 



yr T{2k-t,)T{k + S, + l) , ^,+2^,,+,, 

^^^^^T{k - tj)T{2k + s, + 1) 



(31) 



G(fc+1)^ 
G(2fc + 1) 



E 



rank A 



1} 



r(2fc-tj) 



^(2fe^ 



-1) 



s,!ij!(l + .si + ij) 



.1^1 



dxd 



have equal coefficients of for r < 2A; + 1. For a fixed r, the coefficients of for 
low values of k can be meromorphically continued into each other. 

Furthermore, by using Cauchy's Lemma, one can switch to an expression involv- 
ing products instead of determinants (i.e. a hypergeometric expression). 

Proof. For Expression (29), we proceed essentially by substitution into Equation (22), 
and looking at terms of order N'' A gain, Cauchy's Lemma is used repeatedly 
to reorganize determinants. 

For Expressions (30) or (31), we reorganized yet again the determinants using 
Cauchy's Lemma into a form corresponding to Formula (16). We also summed 
over partitions A instead of summing first over their rank d then their Frobenius 
coordinates ~s , t . 

For a fixed r, both sides indeed admit meromorphic continuations in fc, which 
are equal by Carlson's Theorem [AAR99, Theorem 2.8.1, p. 110]*^. Indeed, the LHS 
is shown to admit a meromorphic continuation in k using a Pochhammer contour. 
The RHS' meromorphic continuation is already written in Expression (30), if we 
admit that what is meant there is the value of the meromorphic continuation in k 
evaluated at k. The difference of the two sides satisfies the hypotheses in Carlson's 
Theorem, in that its value is at integers, it is of exponential type, and type < tt 
along axes parallel to the imaginary axis. The author has shown similar statements 
in his thesis. 

It is probably good to insist that the meromorphic continuation of r|fc^t*)py2fc+s'+ij 
to the left has to be taken very carefully and cannot be obtained by just plugging 
in values of k, once k < t. We will discuss similar issues later, in Section 9. □ 

We now aim to replace the determinant left in Equation (31) by a friendlier 
expression, a rational function of k. 

6. General shape of {M){2k,r), \M\{2k,2h) and \V\{2k,2h) 
We now prove Theorem 1. 



^See also the author's thesis [Dch06b]. 
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Proof. By Equation (29) , we know that (for fixed r and as meromorphic functions 
of A:) 

G(fc+1)2 



-{M){2k,r) 
r\ 



G{2k + 1) 



J2 C{d,^,t)\H 



\dxd ' 



l<d 



with C{d,~s, t ) G Q, while for s and t fixed (and non- negative, of course), Equa- 
tion (27) indicates that H^''^^^ is a rational function of k: 



(32) 



-^k,s,t 



l + s + t 2k 

i— — t 



n 



This already shows that we have a rational function of k and that the degree of its 
numerator equals the degree of its denominator. Equations (32) and (29) together, 
along with the fact that H^''^^^ = f{-k,t,s j-g^ consequence of Equation (27), explain 
why Xr is even. 

In order to determine the YrS a bit better, we need to investigate possible de- 



nominators in the terms of \H 



Idxd' 



If a is positive, Lff 



Idxd 



will have a 



factor of (2k + a)""^''^ in its denominator if and only if a is odd (because there is 



cancellation in Formula (32)) and all of si, 
happen, we need 



(33) 



eta (^^) 

> aa{r)+ ^ (a + i 



, Sa^ are greater than a. For this to 



1)+ £(*-!), 

i=l 



where the inequality is obtained by taking as small as possible values for d (i.e. 
aa{r)), for the Si's (while requiring them to be different and greater or equal to a) 
and for the ti's (all different). 

We turn this inequality around and get 



aa{r) < 



—a + + 4r 



The case of a negative is the same, exchanging the roles played by s and t . 
Finally, the constant D{r) ensuring that both and Yr are monic can be found, 
thanks to Equations (29) and (32), taking limfe^oo^ 



(34) D{r) 



E 



i<d 

d+'Z(si+ti)=r 



1 



1 



dxd 



1 



where this last equality is left to the reader. 

Actually, this last equality is enough to also guarantee that Xr{u) and Yr{u) 
both have integer coefficients: just substitute for i/*^'*'* in Equation (29) 



1 

1 + S + t 



n 



1 

2k + i 



s+t 



k^h.k^ 



for the appropriate (integer) /i^s (in particular, hs+t 
This proves Equation (2), at least for large k. 

Meromorphic continuation has already been obtained in Theorem 15. 



□ 
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Theorem 16. For ft, G N, there are polynomials X2h,X2h, with integer coefficients 
and 6BgX2h = degX2/i > degX2h such that as meromorphic functions of k, 

where Y^^u) is as defined in Theorem 1. 

Moreover (hut this is conjectural) , the numerators are additionally monic poly- 
nomials^ when C{h) = and C{h) ~ , md degX2;i — degX2h — 2/i. 

Proof. For fixed integer r and large integer k, most of this follows immediately from 
Equations (7) and (9), combined with Theorem 1. 

The fact that deg X2h < deg X2h for instance is a consequence of 

(M)(2.,.-,^, f_;\'G(.+ l)= 
which we use in Equation (10): 



2 J G(2k + l) 



2'i /r,i.\ / ■ \ j / ■ \ 2h-j 



E 

J=0 



We can similarly show that if it exists, C{h) = The constant C{h) is more 
mysterious, and involves lower order terms in k of Equation (32). 

The meromorphic continuation is obtained as in the proof of Theorem 1. □ 

Remark. Unfortunately, within their degree restrictions, the Xr{u), X2h{u) 
and X2hiu) polynomials still look utterly random. We merely have an expression 
for them as a sum of determinants of rank d < ^/r (resp. 2ft,). This expression is 
relatively quick and allows at least to compute a few of those polynomials. 



7. Computational data 

7.1. The polynomials Xr{u), X2h{u) and X2h{u)- We present our data for 
{M){2k,r) first, in Table 2, followed by data on \M\{2k,2h) in Table 3 and fi- 
nally on |V|(2fc,2ft) in Table 4. Everything extends numerical results previously 
published, for instance in [Hal04, Hal02a] (but those rely on [Hug05]) or [CRS06] 
(which is limited to A; — ft). Extended versions of those tables are also made 
available in the source of this arXiv submission or (possibly more) at [Deh07]. 

To obtain those tables, we have implemented Equation (30), which is the most 
computationally accessible version of the formulas available in Theorem 15. A 
Magma implementation of this algorithm is also part of this arXiv submission. 

7.2. The roots of X2h{u). It has been suggested before, based on limited numeri- 
cal data, that the polynomials X2h{u) have only real roots. In fact we list in Table 5 
the number of real roots and degree for each such polynomial. One quickly observes 
that X42{u) (of course!) is actually the first polynomial to break the initial fluke 
and have non-real roots. It is not clear at this point if this is related to a similar 
observation on the last line of [Hal02a] and throughout [Hal04]. 



This is the normahzation we will keep later, when discussing data about those polynomials. 
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Table 5. The degree and the number of real roots of X2h- The 
hs for which there are non-real roots are highlighted. 



h 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 




deg{X2h) 








2 


2 


4 


6 


8 


8 


12 


14 




# real roots 








2 


2 


4 


6 


8 


8 


12 


14 






h 


11 


12 


13 


14 


15 


16 


17 


18 


19 


20 




deg{X2h) 


16 


18 


20 


22 


28 


28 


30 


34 


36 


38 




# real roots 


16 


18 


20 


22 


28 


28 


30 


34 


36 


38 






h 


21 


22 


23 


24 


25 


26 


27 


28 


29 


30 




deg{X2h) 


44 


46 


48 


50 


54 


56 


62 


64 


66 


72 




# real roots 


40 


46 


44 


46 


54 


52 


58 


60 


62 


68 





-^42 (u) = " - 12302ti + 69239935U - 236610412148u 

+ 549459541784707u^'^ - 919748248913270486u^* + 1148989069656897835213ii^^ 

- 1094474723973849448826480u^° + 805533314533281755701371226ii^^ 

- 46154192896771811025394423705211^'' + 206514429127544387915748094513446u^^ 

- 7211944111833986997212154158707692011^^ 

+ 1957719645769350260302671962483440450211^" 

- 4099121776759328236737053383626986012604u'^*' 

+ 65417072796093709686120314825046272081985011^'' 

- 7821250373476711537975831731977492624380017611^* 

+ 683698000800342857229690081485643432100615518911^^ 

- 42202825088622350114236559209834534385071085746211^° 
+ 17476800084974190439148752639441918166326024419531u^ 

- 448540393629268182677088044978029477583305447285620ti'' 
+ 625352693721064232359698456539459340167253970973077511* 

- 3701308775622899343826682746064337776289455085124875011^ 
+ 36216052456609571501642100973941635690472733838765625. 

This polynomial has four non-real roots (±18.8631835 ± 0.0090603i) that show 
up at once, since they would have to come in pairs of conjugate pairs by evenness 
of X2h{y')- One could wonder why non-real roots show up so late, and if there is 
actually a good reason for this. 

Fact. The polynomials X2h{u) tend to have many, but not all, of their roots real. 
For instance, for high h, X2hiu) has one root very close by to every odd integer 
between h and 2h. 

We first present graphical clues for this fact in Figure 1 , which depicts the posi- 
tion of the real roots for ft, = 1 to ft. = 30. It thus omits the complex roots. 

We now explain the fact. It helps at this point to remember that X2h{u) is 
obtained by summing various Xr{u) for r < 2ft, which are themselves obtained 
from Equation (29), for instance. Furthermore, the summand in that Equation 
associated to d,s,i (with r = d+J2i Si + ti) will have poles (as a function of u = 2k) 
at the odd integers a such that — si < a < ii (this uses Lemma 5 to expand the 
determinant in H^-^^'^^s). For each pole a, there are a few summand where this pole 
comes with multiplicity exactly aa{r), but for most others the multiplicity is lower 
(see Equation (33)). So if we sum all of those terms, and multiply by Y2h{u) (the 
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Table 6. The largest root of X2h{u) 



h 


largest root of X2h{u) 


difference with 2h — 1 


log. difference 


1 




no root 


no root 


no root 


2 




no root 


no root 


no root 


3 


3.0000000000000000000 


2.000000000 


0.69315 


4 


5.7445626465380286598 


1.255437354 


0.22748 


5 


8.2448923938491831987 


0.7551076062 


-0.28090 


6 


10 


568920444013080343 


0.4310795560 


-0.84146 


7 


12 


769459455674733521 


0.2305405443 


-1.4673 


8 


14 


886048429155973920 


0.1139515708 


-2.1720 


9 


16 


948550444560344620 


0.05144955544 


-2.9672 


10 


18 


978943770872905688 


0.02105622913 


-3.8606 


11 


20 


992206162055831068 


0.007793837944 


-4.8544 


12 


22 


997383184072186530 


0.002616815928 


-5.9458 


13 


24 


999198051064882757 


0.0008019489351 


-7.1285 


14 


26 


999774030173017860 


0.0002259698270 


-8.3951 


15 


28 


999941044846106152 


5.895515389 x 10"^ 


-9.7388 


16 


30 


999985671005722891 


1.432899428 x lO"'^ 


-11.153 


17 


32 


999996738730003824 


3.261269996 x 10"*^ 


-12.633 


18 


34 


999999301847217917 


6.981527821 x 10^'^ 


-14.175 


1 Q 


36 


999999858891343014 


1/111 nsRK7n \r 1 n~''' 


1 K 77 A 
-ID. 1 / 4 


20 


38 


999999972983353984 


2.701664602 x 10~^ 


-17.427 


21 


40 


999999995085836086 


4.914163914 x 10^^ 


-19.131 


22 


42 


999999999148595422 


8.514045781 x 10^^" 


-20.884 


23 


44 


999999999859167358 


1.408326421 x lO^i" 


-22.683 


24 


46 


999999999977712180 


2.228782021 x 10^" 


-24.527 


25 


48 


999999999996618870 


3.381129731 x lO'^^ 


-26.413 


26 


50 


999999999999507453 


4.925468142 x 10"" 


-28.339 


27 


52 


999999999999930988 


6.901186254 x 10"!^ 


-30.304 


28 


54 


999999999999990686 


9.313971788 x lO'^^ 


-32.307 


29 


56 


999999999999998787 


1.212486889 x lO-^^ 


-34.346 


30 


58 


999999999999999847 


1.524414999 x lO^^^ 


-36.420 



common denominator) to obtain X2hiu), a vast majority of terms factor a (m — a) 
out. We thus have an expression of the form 

X2h{u) - {u-a)Pi{u) + P2{u), 

where the coefficients of Pi (u) are expected to be much bigger than the coefficients 
of P2(m) (simply because much more terms are summed to obtain Pi{u) than P2(w)). 
Hence, we should expect X2h{u) to change sign when u travels along the real axis 
from below a to above a (because |Pi(a)| >|P2(a)| and (u — a) changes sign) and we 
know that a root will be around u = a. This is especially true if a > r/2, because 
the restrictions impose then si > a > S2, and as a consequence ar{a) — 1 and 
the phenomenon described just now is accentuated. We present in Table 6 some 
numerical data associated to this phenomenon. 

It is obvious from Figure 1 that a lot is yet to be understood about the polyno- 
mials X2hiu). For instance, it is not cleax if asymptotically in h there is a positive 
proportion of real roots. 
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8. Alternative expressions 



8.1. Using Macdonald's ninth variation of Schur functions. Define, as in 
[NNSYOl] and [Nou04], and similarly to [Mac92], 



(35) 
with 



l{X)xl(X) 

■ {R + k-iy.k\' 



We first prove that this variation of Schur functions satisfies a Giambelli identity. 

Proposition 17. Let X be a partition and Frohenius coordinates, of rank d. 

Then, 



dxd 



r(_R+s,+l) 



dxd 



Note how this provides a second determinantal expression for this variation of 
Schur functions, but with a matrix of different rank. 

Proof. We intend to use Exercise 3.21 in Macdonald's book, but to show that the 
exercise applies, we need to prove: 



det 



:iR) 

1 




'p+2 



1 



■(R-2) 7{R~3) 



V 
1 







liR-q) 

"'p+q+ 

t{R-q) 



1 



{q+l)x{q+l) 



T{R - q) 



p\q\{l+p + q)T{R + p+l)' 
This can be shown by expanding the determinant along the last column to obtain 



^ip\q) 



Subtract the LHS from the RHS, proceed by induction on q, factor out p^^^p^\) and 
the result then follows from the following equalities, for p and q positive integers: 



(-1)« 



E: 



(-1)' 



(p+q+iy. f-^ ilpl{q ~ iy.{p + q - i + I) p\q\{l+p+q) 

(-1)« , P + q + q 2Fi(i-?_pip;l) 1 



IV. 



(p + q)ip + q + ^)pW- 



p\q\{l+p + q) 



-p-q. I) 
p-q ' ' 



= 0, 



{p+q+l)\ {p + q){p + q+l)p\q\ 

{-lY I q (g-i)! 

{p+q+l)\ {p + q){p + q + l)p\q\ (1 - P - 

the last one being a consequence of Gauss's Hypergeometric theorem. 

The theorem now results directly from Exercise 3.21 in Macdonald's book. □ 
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In essence, this Proposition allows us to switch from a Giambelli-type expres- 
sion to a Jacobi-Trudi expression. It immediately leads to a simplified version of 
Theorem 15. 

Theorem 18. With G(-) the Barnes G-function, and sx defined as in Equa- 
tion (35), 

in the sense that their coefficients of are equal for fixed r and large enough k so 
the coefficient in the LHS is defined. 

8.2. Imitating the Cauchy identity. We can also give an alternative for the 
expression in Equation (31), proceeding as in Gessel's theorem in its lead up to the 
Cauchy identity (see [TWOl]). This uses Theorem 18. 

Theorem 19. 

G(fc + 1)2 v^,,^,,^,^ ^^{yzf 



iim^detl5:;.,_.([i'^])er"+^)z'-^- 



r>0 

= lim det 

n — >oo 

= lim det 



lim det > 

fi — » on I ' ^ 



l>0 



l^i + k~l\ (2fc-n + j-l)! 

fc-1 ) {l-jy.{2k-n + l-iy.^ 



in the sense that their coefficients of z^ are equal for fixed r and large enough k so 
the coefficient in the LHS is defined. The factorials on the last line should really be 
evaluated in groups, to give if I < j, and p ( 2k- +"k j ) ! otherwise. 

Note that this can be truncated significantly when we are after only 



(M)(2fc,r)M: 

0<r<S 



for a finite S (i.e. when we are computing the head of the sequence of polynomials): 
we can drop the limit in n and settle for a sufficiently big n instead, and then cut 
the matrices in their infinite directions as well. 

In Gessel's Theorem, in order to get to the other side of the Cauchy identity, 
one would then observe that the matrix on the last line is Toeplitz, and then use 
Szego's theorem. Of course, that fails here because the matrix on the last line is 
not Toeplitz. 

9. The result of Conrey and Ghosh 

As explained in the introduction, Conrey and Ghosh's theorem [CG89] that 
J{2, 1) = immediately leads to a conjecture that |V|(2, 1) = as well. Our 
main concern is that we only know | V|(2fc, 2/i), for integer h, through Equations (7) 
and (9) (while we would need h = 1/2). 

We offer in Figure 2 one way to circumvent this problem. The idea is to com- 
pute for each fixed integer h the values of the mcromorphic continuation in k of 
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Figure 2. The real part of the situation in the Conrey-Ghosh 
case. The circle at (1, 1/2) indicates the point for which the value 
of J {2k, 2h) is coveted. The dots indicate the locations where Ex- 
pression (30) applies, and the crosses indicate the points to which 
that expression is meromorphically continued (for a fixed h, i.e. 
horizontally) thanks to Expression (36). Note that for fixed inte- 
ger h, this continuation hits a pole when crossing the dashed line 
(and many more before reaching k — 1, as /i increases: see Fig- 
ure 1). 



{M){2k, 2h) at fc = 1 (i.e. at the crosses). This should be enough to know through 
Equation (7) any value of the form \ \ (2, 2h), which could then finally be used to 
meromorphically continue |V| (2,2ft,) to h ~ 1/2. 

Getting the meromorphic continuation of Equation (30) to fc = 1 is quite subtle. 



Proposition 20. Define {A4) {2, r) as the meromorphic continuation in k of (A4){2k,r), 
evaluated at k = \. Then, the exponential generating series of (Al)(2,r) is given 
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hy 

(36) ^(Al)(2,r)M: 



r>0 



^ ^ s,!i,!(i+s,+t,) -In 



where i^(0) = w/ien t = Q, v{t) = 1 w/ien t>2. The value v{\-) is free to choose. 

Proof. When looking for the analytic continuation in k, most of the formulas we 
have found so far are misleading. For instance, in light of the remark in footnote 5, 
one could think that the sums over partitions A in Expression (31) or Theorem 18 
immediately reduce when A: = 1 to sums over partitions A of length 1, i.e. partitions 
indexed by a single variable. However, in those cases, the other factor in the 

(2k) 

summands (i.e. for instance s\ in Theorem 18) might actually be undefined if 
we take A: = 1 (in that particular case, when 1{X) > 3^°). 

We can get a better intuition through Expression (30) , which we use as a basis of 
our proof. We are clearly required to find the meromorphic continuation to fc = 1 
and for fixed s, t > G N of 

r(2fc-i) r(fc + s + i) 

T{k-t) ' T{2k + s+l)' 
The second factor is certainly not a problem and immediately gives -^rp2- ^^^^ 
first factor, we have to look at limfe^i ^-^"^^Zq for t > 0. Pick any integer a such 
that I + a — t > 0. Then, using the functional equation for T, we have 

lim ^''^'^ ~ - lim + {k-t){k-t + l)---{k-l)---{k + a~t) 

fe™ T{k~t) ~ fc™ r{k + a-t+l) ' (2fc - t) • • • (2fc - 2) • • • (2/s + a - i - 1) 

Note that the terms 2k^2 ^^^y appear if t > 2. In that case we get 

,. r(2fc-<) , r(2fc + a-i) fc-1 

lim / = (1 - t hm — ^ '-- = - 1 - t 

k^iT{k-t) ^ ' k^iT{k + a-t+l) 2k-2 2^ 

and in the case t < 2 the factor of 2 is missing. □ 
One can also check that the values recovered using Proposition 20 agree with 

JC (2) 

the values obtained using Yij2j ^^'^ thus Theorem 1. 

For completeness, we give the beginning of the sequence of Xr{2)s, for r = 1 
to 15: 

1, 2, 0, 18, 50, -6540, -11760, 852180, 1228500, 590126040, 558613440, 

- 39273224760,455842787400,5775116644337040, 14904865051876800 

Unfortunately, we fall short of actually finding the full meromorphic continuation 
of {A4){2,r) and have to leave this for a further paper. 

10. Conclusion 

The initial goal was to compute the (TW) (2fc, r), \M\{2k,2h) and \V\{2k,2h) 
more effectively than previously done. 

We feel that we have achieved this goal, since we have been able to shed some 
light (for instance in Theorem 1) on the structure of the results. This structure 
(rational functions with known denominators) underlines tables already available 
in [Hug05] or [CRS06]. We have also been able to use these results to obtain better 



^"or equivalently when ti > 2 if A = 
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algorithms to compute those rational functions, thereby extending the data that was 
available. Much of that data is now available in the source of the a rXiv submission, 
or at [Deh07]. As a corollary we have shown that for large(r) h the roots (in k) of 
I V|(2fc, 2h) cease to all be real, a fluke only for the small-/i cases available previously. 

However, we have not obtained a formula for all |V|(2fc,r). In particular, we 
cannot recover the value of |V|(2, 1), which can be conjectured from Conrey and 
Ghosh's result for J (2,1). 

Those methods should also give more general moments, for instance for expres- 
sions of the form 



\Zu{9i)\ 



2k 



Z'u{()2) 



Zu{02) 



or 



(37) 



\Zu{ei)\ 



2k 



Z'{j{ 



Zu{&2) 



U(JV) 



U(N) 



An expression for those two extensions in the shape of Equation (15) would defi- 
nitely be available (for instance, in the case of Expression (37), we would most likely 
have to compute the equivalent of Equation (15) by summing over /x e (2N+)'"). 
However, the second part of the computation, the part covered here by Proposi- 
tion 10, would probably be significantly worsened. 
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